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SECTION  I 
INTRODUCTION 


In  many  applications  it  is  necessary  to  place  a  protective  cover 
over  an  antenna.  Such  a  cover  is  commonly  referred  to  as  a  radome  and 
the  conventional  approach  has  been  to  use  solid  or  laminated  dielectric 
materials  for  these  covers.  Recently,  substantial  effort  has  been 
directed  to  a  study  of  metallic  radomes. 

Metallic  radomes  possess  many  inherent  advantages  over  conven¬ 
tional  dielectric  radomes.  Some  of  these  are 

1)  Elimination  of  the  precipitation  static  (p-static)  noise, 
which  may  cause  certain  radars  and  electronic  equipment 
to  malfunction, 

2)  Inherent  lightning  protection, 

3)  Reflection  from  a  thermonuclear  flash, 

4)  Better  shielding  against  low  frequency  EM-pulses  since  a 
metallic  radome  in  general  can  be  regarded  as  a  bandpass 
filter  with  high  attenuation  for  low  frequencies, 

5)  Potentially  better  laser  hardening,  and 

6)  Potentially  higher  mechanical  strength. 

Periodic  surfaces  are  employed  in  the  design  of  metallic  radomes. 
Previous  designs  of  periodic  surfaces  have  been  used  for  band  rejection 
of  an  incident  signal  by  utilizing  dipole  arrays  and  they  have  been 
used  as  narrow  bandpass  filters  by  utilizing  slot  arrays.  The  objec¬ 
tive  of  this  study  is,  through  the  use  of  planar  periodic  slot  arrays, 
to  design  a  passband  filter  with  the  largest  bandwidth  to  date.  The 
bandwidth  is  considered  to  be  that  range  of  frequencies  in  which  the 
field  transmitted  through  the  structure  has  a  loss  of  less  than  1  dB. 

In  some  applications  a  typical  constraint  is  given  as  3  dB  but  this 
would  mean  50%  of  the  energy  is  lost.  This  lost  energy  is  not  neces¬ 
sarily  a.teorbed  but  can  be  reflected  in  directions  different  from  that 
on  the  m®n  beam.  This  could  have  the  undesirable  effect  of  increasing 
the  sidejfobe  level  of  the  antenna. 


The  type  of  metallic  radome  to  achieve  the  objective  is  found  to 
be  an  extension  of  previous  results.  From  the  results  in  [1],  it  is 
determined  that  the  type  of  element  in  the  planar  periodic  slot  array 
be  a  generalized  three-legged  element.  This  type  of  element  allows  for 
a  closer  packing  of  the  elements  and  also  allows  the  elements  to  be  un¬ 
loaded.  From  the  results  in  [2],  it  is  determined  that  a  skewed  grid 
(i.e.,  interlace  design)  should  improve  the  bandwidth.  From  the  results 
in  [3],  it  is  determined  that  two  planar  periodic  slot  arrays  placed 
in  parallel  results  in  an  increase  in  bandwidth  with  more  control  of 
the  shape  of  the  transmission  curve.  This  results  in  a  flat  response 
in  the  pass-band  region  and  greater  attenuation  outside  the  pass-band 
region. 

There  must  be  a  dielectric  structure  to  be  used  to  separate  and 
to  support  the  two  arrays.  From  the  results  in  [4]  it  is  determined 
that  dielectric  layers  are  needed.  These  layers  are  used  to  provide 
a  constant  bandwidth  for  the  varying  angles  of  incidence.  From  the 
results  in  [5],  it  is  determined  that  three  dielectric  layers  should 
be  used.  The  dielectric  layers  serve  the  dual  function  of  physical 
support  as  well  as  stabilizing  the  bandwidth.  The  outer  dielectric 
layers  produce  a  constant  bandwidth  and  the  middle  dielectric  layer 
provides  the  proper  coupling  between  the  arrays. 

All  of  this  leads  to  the  structure  of  Figures  1  and  2,  that  being 
a  biplanar  slot  array  of  generalized  three-legged  elements  imbedded  in 
three  dielectric  slabs.  The  slabs  have  relative  permittivity,  erm,  and 
thickness,  dm,  where  m  refers  to  the  dielectric  media.  The  theory  is 
described  in  the  next  chapter  and  the  actual  transmission  curves  and 
input  variables  are  determined  in  Chapter  III.  The  conclusions  are 
presented  in  Chapter  IV. 

The  approach  used  in  this  report  is  to  use  the  previous  solutions 
as  a  guide  and  to  apply  them  with  a  view  toward  obtaining  as  large  a 
bandwidth  as  possible.  Since  much  of  the  theory  is  contained  in  detail 
in  these  prior  studies  by  other  authors,  it  is  not  repeated  here.  Only 
the  final  equations  are  presented  so  that  the  casual  reader  can  observe 
what  is  being  achieved.  This  also  allows  the  achievements  to  be  pre¬ 
sented  with  a  minimum  of  clutter.  The  serious  reader  will,  of  course, 
refer  to  these  original  and  rather  elaborate  reports. 
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Figure  2.  Three-legged  slot  array  geometry. 
L.  A.  denotes  Leg  Angle. 
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StCTION  /I 

DEVELOPMENT  OF  THE  THEORY 


Generalized  three  legged  elements  vonsist  of  three  monopoles  con¬ 
nected  together  at  a  single  point  each  of  arbitrary  length  and  direc¬ 
tion.  The  lengths  being  ft1*1,  ft2*1,  s.3*1  and  the  directions  being  unit 
vectors  p1 »i ,  02, i,  p3,i  where  the  superscripts  refer  to  the  leg  number 
and  array  index,  respectively*  The  parameters  of  the  generalized  three 
legged  element  are  illustrated  in  Figure  3.  The  lengths  and  unit  vec¬ 
tors  will  remain  fixed  in  each  array  but  mu*/  differ  between  arrays  for 
the  remainder  of  this  report.  It  has  been  shown  in  earlier  work[6] 
that  the  voltage  distribution  along  the  reference  slot  in  array  i  is 
the  sum  of  two  modes: 

1)  the  symmetric  mode,  Vs*1 (t)  i  -  1 ,2 

2)  the  asymmetric  mode  Va » 1  (jt ) 
as  shown  in  Figure  4. 


Figure  3.  Generalized  three  legged  element 
showing  critical  parameters. 


*  Note:  The  superscript  is  sometimes  enclosed  with  a  parenthesis. so 
that  itjs  not  confused  with  the  power  of  a  variable  (i.e.,  p^21 
is  not  p  squared). 
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Figure  4.  The  voltage  modes  on  a  generalized 
three  legged  element. 


To  determine  the  transmission  through  this  structure,  the  proce¬ 
dure  is  similar  to  that  of  an  earlier  report[7].  This  consists  of  the 
following  three  parts: 

A.  determination  of  the  symmetric  current, Is,in,  and  the 
asymmetric  current,  Ia'in  induced  by  the  incident  H-field, 

B.  determination  of  the  symmetric  and  asymmetric  voltage 

modes,  Vs>2(t)  and  respectively,  where  the  super¬ 

script  2  refers  to  the  second  array, 

C.  determination  of  the  transmitted  H-field  reradiated  by 
the  two  voltage  modes  above  of  the  second  array. 


A.  Determination  of  the  induced  currents  is,in  and  la,in 

Let  the  configuration  shown  in  Figure  1  be  exposed  to  the  inci¬ 
dent  plane  wave  whose  magnetic  field  is  given  by 


H]nc(R)  =  HjnC 


-jeiR-S] 

e 


where 

Si  is  the  direction  of  the  incident  Diane  wave  signal 
in  medium  1  (as  shown  in  Figure  5), 
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Figure  5.  Coordinate  system  for  incident  fields. 


f$l  is  the  propagation  constant  of  medium  1, 

R  is  the  position  vector  for  the  point  of  observation 
for  the  incident  field  (as  shown  in  Figure  6). 

A  slot  can  be  considered  as  a  magnetic  element  mounted  directly 
in  front  of  an  electrically  perfectly-conducting  ground  plane[8].  Be¬ 
cause  of  this  electric  screen  effect,  the  incident  plane  wave  will. only 
induce  ls»in  and  Ia>1n  in  array  1.  Also  the  induced  currents,  Is»in 
and  Ia,1n,  are  independent  of  whatever  exists  behind  array  1.*  Thus 


*  The  effect  of  the  structure  behind  array  1  is  contained  in  the  mutual 
coupling  terms  between  the  two  arrays. 


ybm 


OBSERVATION  POINT 
OF  INCIDENT  FIELD 

Figure  6.  Structure  used  to  define  certain  variables 
needed  to  determine  the  admittances  in  slab 
m  bounded  by  the  planes  y=bm  and  y=bm+^. 

the  induced  currents  are  of  the  same  form  as  Equations  (11)  and  (12) 
in  [9).  The  symmetric  and  asymmetric  current  are: 


Is’1n  =  [H]nc(R)  •  in1  P*U  aT2/1(0,0)  +  Hj^(R)  •  ^  / 

-je?(R(1)-R)*s2 

11^2/1  3  e  0  ) 


,slt 


ja.iii  =  [jrinc(R)  .  ~  palt  t  (0,0)  +  Hinc(R)  ■  n,  PaU 

I  *  I  1  1  2/i  1  m  1  li 


-jMR^-RKo 

h^2/1^*®^  e  (2) 
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where 


in1  is  a  unit  vector  orthogonal  to  the  plane  of  incidence 
in  medium  1  (see  Appendix  C), 


„n^  is  a  unit  vector  parallel  to  the  plane  of  incidence 
and  orthogonal  to  the  direction  of  propagation  in 
medium  1  (see  Appendix  C), 

is  the  position  vector  for  the  reference  point  of  the 
reference  element  in  array  1  (see  Figure  6)» 


{»> 


is  the  composite  symmetric  transmitting  pattern 

factor  for  array  1  in  medium  1,  orthogonal  and  parallel 
components. 


^ijP]  is  the  composite  asymmetric  transmitting  pattern  factor 
11  for  array  1  in  medium  1,  orthogonal  and  parallel 
components , 


lii' 


1^^ (0,0)  is  the  transformation  function  for  dielectric 
media  2  normalized  to  media  1.  (normalized  T-factor). 


The  pattern  factors  are  the  far  field  patterns  due  to  a  voltage 
distribution  on  an  element  of  the  array.  They  represent  the  relative 
magnitude  of  the  plane  waves  propagating  in  certain  directions.  See 
section  (J2,  The  Pattern  Factors,  for  defining  equations. 

The  transformation  function  represents  the  transformation  of  the 
field  as  a  result  of  the  dielectric  interfaces.  The  plane  waves  will 
be  partly  reflected  and  partly  transmitted  at  each  interface.  The 
transformation  function  sums  up  this  effect  of  multiple  reflections  of 
the  waves  and  transforms  tne  field  from  one  dielectric  media  to  another 
see  section  B3,  The  Transformation  Functions,  for  defining  equations. 


B.  Determination  of  Vs,^(t)  and  Va,^(si.) 

Certain  quantities  needed  for  the  evaluation  of  Vs»^(«.)  and  Va,2U) 
are  evaluated  in  the  following  sections. 

Bl.  Determination  of  Admittances 

As  developed  in  [10],  the  impedances  for  dipole  arrays  was  defined 

as 


zi',i  E  _  Voltage  induced  in  the  single  test  element  i1 _ 

Terminal  current  of  the  reference  element  of  the  array  i* 


♦This  does  not  imply  that  the  currents  in  the  remaining  slots  are  zero. 
In  fact  they  are  related  to  the  reference  element  by  Floquet's  theorem. 
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1 


Through  use  of  duality  (i.e.,  the  induced  voltage  becomes  the  induced 
current  and  the  terminal  current  becomes  the  terminal  voltage),  the 
admittance  for  slot  arrays  is  defined  as 

Yi',i  .  _  Current  induced  in  the  single  test  element  i1 _ 

Terminal  voltage  of  the  reference  element  of  the  array  i 

Now  each  voltage  mode,  Vsl(t)  and  Vai(t)  of  each  array  induces  a  sym¬ 
metric  current  in  the  reference  element  of  each  array.  Hence  there  is 
a  total  of  sixteen  admittances. 

Using  the  dual  of  Equation  (44)  in  [11]  with  some  change  of  nota¬ 
tion  and  generalizing  to  the  medium  m  results  in 

-j6m(R(i')-R(i))-r 
e  m  m 

r 

my 


yAi'.Bi  _ 


2DXDZ 


I  I 

k=-<x>  n= 


r  pAi '  t  pBi  T  / .,  ( i )  d  _|i(  i 1 )  \  +  pAi't:  pBi  T  fu(i)  d  _u(i'))-| 
LiHm  iimvu  ,am  u  '  ir m  «Km  h 'm'0  ,am  u 


(3) 


where 


A,B  are  dummy  superscripts  that  refer  to  the  symmetric  or 
asymmetric  composite  pattern  factors, 

.x-P™1  *V-i^mlare  defined  in  the  section  on  Pattern  Factors, 

Ym  is  the  characteristic  admittance  of  dielectric  media  m, 

DX,DZ  are  the  interelementspacings  between  adjacent  slots 
measured  in  the  x  and  z  direction  respectively, 

Bjj,  is  the  propagation  constant  of  media  m, 

u(i),u(1")  are  defined  in  Figure  6, 

R ( 1 ' )  is  the  position  vector  of  the  single  test  element 
reference  point, 

)  is  the  position  vector  of  the  reference  point  of  the 
reference  element  of  array  i, 

x  Tp/u^J.dpj-u^'))  is  defined  in  section  B3,  Trans- 
^  ir 

formation  Function. 
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The  unit  vector,  rm,  indicates  the  directions  in  which  the  bundle  of 
plane,  inhomogeneous  waves  are  being  scattered  by  the  array.  Using 
Equations  (D2)  and  (D3)  in  [12]  together  with  results  developed  in  [13]  , 

rm  for  a  skewed  array  is  found  to  be 


x 


+ 


nAzXm  \ 
°x  DxDz  ) 


yr, 


my 


+  z 


(4) 


where 


rmy 


+ 


nAzx  \ 

_ 51) 

DXDZ  ) 


+ 


nx. 


2.1/2 


(5) 


and 


az  is  the  interlace  spacing  in  the  z  direction,  and  is  defined 

as  the  distance  between  two  slots  that  are  adjacent  in  the 

x  direction  (see  Figure  2), 

Am  is  the  wavelength  in  media  m. 

It  can  be  shown  that  results  obtained  interlacing  in  the  x  direction  are 
equal  to  those  obtained  from  the  analysis  for  interlacing  in  the  z  direc¬ 
tion.  Also  the  solution  for  geometries  interlaced  in  both  directions  at 
once  can  be  evaluated  by  redefining  Dx  and  Dz  to  new  values. 

Now  r  may  either  be  real  or  imaginary.  When  it  is  real,  r^y 
represents  a  wave  propagating  away  from  the  array,  hence  the  plus  sign 
must  be  used.  When  it  is  imaginary,  r^y  represents  a  wave  which  attenu¬ 
ates  as  it  moves  away  from  the  arrajf,  nence  the  negative  sign  must  be 

used.  For  k=n=0,  the  equation  for  ?m  reduces  to  sm  (and  the  specular 

di recti  on) . 


The  self  admittance  is  closely  approximated  by  the  mutual  coupling 
between  the  single  test  element  located  w/4  (w  is  the  width  of  the  slot) 
away  from  the  reference  element  of  the  array [14].  Assuming  that  the  x 
components  and  the  z  components  of  RU1)  and  Ru)  are  equal. 


)  ‘  rm  =  4  rmy 


(6) 


For  mutual  admittances,  the  single 
element  of  the  other  array.  Using 
assumptions  as  before,  results  in 


test  element  is  assumed  to  be  an 
the  given  structure  with  the  same 
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d3r3y 


(7) 


Using  Equations  (6)  and  (7)  will  assist  in  determination  of  the  self 
and  mutual  admittances.  Appendix  A  contains  a  complete  list  of  all 
the  admittances  for  the  structure  in  Figure  1. 


B2.  The  Pattern  Factors 


The  transmitting  and  non-transmitting  pattern  factor  for  each  leg 
is  found  through  the  application  of  duality  to  Equations  (D14)  and  (40) 
in  [15]  to  give 


pvU(g) 


_J _ 

VvU(0) 


(8) 


and 


Pvi(g) 


_ 1 _ 

Vvi(0) 


VvV)  e' 


dJ i 


g  =  1,2,3 


(9) 


respectively,  where  g  is  the  leg  number  and  i  again  refers  to  the  array 
i ndex . 


There  is  a  symmetric  pattern  and  an  asymmetric  pattern  correspond¬ 
ing  to  the  symmetric  and  asymmetric  voltage  modes  respectively.  Hence 
the  dummy  superscript  1 v 1  for  variable  becomes  an  's'  for  the  symmetric 
voltage  mode  and  becomes  an  'a1-  for  the  asymmetric  voltage  mode. 

Note  that  it  is  not  necessary  to  specify  the  media  for  8  and  r  since 

it  was  assumed  that  the  arrays  are  planar  (i.e.,  p(?)  contains  no  y  com¬ 
ponent).  Using  Equations  (Bl)  and  (B2)  in  [16]  it  can  be  shown  that 

p^^*fm  =  Bq  p^^-rq  (for  planar  elements  only)  (10) 

where  m  and  q  refer  to  the  dielectric  media.  Hence  the  media  subscripts, 
m  and  q,  can  be  eliminated. 

The  form  of  Vvi*U)  and  Vvi(t)  is  usually  assumed  to  be[17] 

VvU(*)  =  VvU(0)  sin8d(^)i-^9'i)  (11) 
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and 


Vvi(0  =  Vvi  (0)  sinf^U9*1-*9,1)*  (12) 

where 

is  the  effective  dielectric  propagation  constant 
(Equation  (C-2)  in  [18]), 

Jief  is  the  effective  length  of  the  leg  due  to  the 
inductance  at  the  ends  of  each  leg. 

This  inductance  effectively  increases  the  length  of  the  legs. 
Vvi*(0)  and  Vvl(0)  are  the  magnitudes  of  the  voltages  for  each  mode 
at  the  terminals  (i.e.,  point  where  legs  join) .  It  must  be  noted  that 
only  the  form  of  the  voltages  Vvi*()i)  and  Vv1U)  is  used  to  calculate 
the  patterns  since  the  magnitude  at  the  terminals  divides  out.  The 
magnitude  of  each  voltage  for  the  second  array  will  be  determined 
later.  It  is  needed  to  calculate  the  transmitted  H-field.  (PartC). 

There  is  a  generalized  composite  pattern  factor  for  the  trans¬ 
mitting  and  non-transmitting  (scattering)  case  for  each  voltage  mode. 

By  inspection  of  Figure  7b  and  the  use  of  Equation  (42)  in  [19]  the 
generalized  composite  pattern  factors  for  the  generalized  three-leg¬ 
ged  elements  are 


P 

(4)  m 


sit  =  (2-(l)ipsU_-(2)ps2t_-(3)ips3t) 


alt  .  <5<2>tp.2*.f(3)ytt)  A 

|J  t|M 


Psi  =  (2p(l)ipsl.p(2)ips2.p(3)ps3)  .  n„ 


Pai  ,  (p(2)pa2_p(3)pa3j  A 
i  ii )  t  ii  I 


*An  alternate  form  for  Vv1(^)  will  be  presented  in  section  III,  Development 
u  Data  and  Results.  The  reasons  for  this  change  are  discussed  in 
that  section. 
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Figure  7.  Showing  excited  modes  in  v i  SCAN  PLANE 

(a)  Loaded  straight  slot  (L.S.S.) 

(b)  Three  legged  direction  vectors 
and  voltages  defined 

(c)  L.S.S.  for  parallel  incident  polarization 

(d)  Three  legged  parallel  incident  polarization 

(e)  L.S.S.  for  orthogonal  incident  polarization 

(f)  Three  legged  orthogonal  incident  polarization. 
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“■  ^~*Tp .  '  "*~fr 


The  factor  of  2  comes  from  the  fact  that  Vs(t)  is  being  counted  twice 
on  leg  a1).  The  negative  signs  appear  because  the  vector  ‘magnetic 
current'  is  going  in  the  opposite  direction  of  the  unit  vectors  of  leg 
directions  as  shown  in  Figure  7b.  All  patterns  have  now  been  specified. 


B3.  The  Transformation  Functions 

The  transformation  functions  (T  factor)  needed  in  admittance  cal¬ 
culations  are  found  from  Equation  (29)  in  [20]  by  changing  from  reflec¬ 
tion  coefficients  for  the  electric  field  to  reflection  coefficients  for 
the  magnetic  field  and  generalizing  the  results  to  layer  m 


,  Tn>{l)’dr,fu(l  ^ 
ill  i 


-j26mu 


(i) 


x 


-j2Bm(dm-u(1,)r, 


my 


y 


“jgmdmrmy 


07) 


This  is  the  generalized  non-normal i zed  T-factor  for  slot  arrays  such 
that  the  array  and  the  single  test  element  are  separate^  by  at  most 
only  one  dielectric  layer  for  u(i')>u(l).  For  u'1'>u(1  )  see  remarks 
at  the  end  of  this  section  For  the  given  structure,  P  .  and 

m  J  iu  m+]  are  the  reflection  coefficients  between  two  dielectric  media 

whose  equations  are  listed  in  Appendix  D.  The  T-factors  used  in  admit¬ 
tance  calculations  can  now  be  computed.  See  Appendix  B. 

It  was  found  that  expressing  the  induced  currents  and  the  trans¬ 
mitted  H-fteld  in  a  normalized  form  results  in  conceptually  simpler 
equations.  The  T-factors  used  in  these  equations  are  hence  normalized. 
Modification  of  Equation  (33)  in  [21]  in  the  same  manner  Equation  (29) 
was  modified  in  [22]  gives 
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(i)Wu(i),dB-u<n)  =  fiV^ 

\  { n  )*  m»l 


(j 


_j2emu^  ^rmy 


)  0 


-j2Bm(dm-u^  ^)r 


-J2®rArmy 


This  is  the  generalized  normalized  T-factor  with  u^1  )>u^)  and  where 
q=m±l  depending  on  whether  calculation  pf  induced  currents  or  trans¬ 
mitted  field  is  involved.  For  u'1/>uvi  )  again  see  the  remarks  at  the 
end  of  this  section. 

The  following  sign  changes  must  be  made  for  the  case  u^W*^  as  is 
shown  in  Appendix  D  in  [23] 


change  from 

subscript  + 
subscript  - 


Determination  of  *5>i  and  V3*1  i=l,2 


subscript  - 
subscript  + 


The  form  of  Vs’1 (t)  and  Va,1(t)  has  been  assumed  previously  (see 
Equation  (12))  in  order  to  calculate  the  patterns.  The  magnitude  of 
Vs’1 (0)  and  Va * 1 ( 0 )  is  now  calculated.  .Denote  the  load  admittances  for 
the  syrmetric  and  asymmetric  mode  by  Yp1  and  Yp1 ,  respectively.  Hence, 
from  Kirchoff's  current  law  (i.e.,  thesum  of  all  currents  entering  a 
node  must  equal  0),  the  load  admittance  multiplied  by  the  unknown  volt¬ 
age  will  equal  the  sum  of  induced  currents. 

Yp1Vs’1(0)  =  Is,in  +  Isls1  +  lslal  +  isls2  +  isla2  (19) 


Ypl va»l (o)  =  ia>in  +  jal si  +  jalal  +  jal s2  +  jala2 


Yp2Vs,2(0)  =  0  +  Is2a2  +  ls2a2  +  ls2s1  +  is2al 
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(22) 


^3 2yd»2^Q^  _  Q  +  jd2s2  ^  |d2d2  ^  jd2sl  ^  jd2dl 


Using  the  equations  for  admittances  developed  in  the  preceding  section 
results  in: 


jS.in 

ja,in 

0 

0 


=  CY] 


where  [Y]  is  defined  as  the  admittance  matrix 


[Y]  = 


Sl+Y[] 

ySlal 

y  S  1  S  2 

yS  1  a2 

al  si 

Ya  1  al  +Yal 

yals2 

yala2 

s2sl 

yS2al 

Ys2s2+Ys2 

ys2a2 

a2sl 

ya2al 

Ya2s2 

ya2a2+Ya2 

All  symmetric  and  asymmetric  voltages  for  each  array  can  now  be  found. 
Since  the  transmitted  H-field  into  the  semi -infinite  space,  y>d4,  is 
desired,  only  the  node  voltages  from  the  second  array,  Vs,z(0)  and 
Va,2(0),  need  be  found.  Using  Cramer's  rule  with  the  determinant  of 
the  admittance  matrix. 


gives 


VS,2(0) 


Sl+Ysl 

ySlal 

jS,in  ySla2 

a  1  s  1 

yd  1  3 1  yd  1 

ja,in  yd 1 a2 

s2sl 

ys2al 

g  yS2a2 

a2sl 

ya2al 

g  Ya2a2+ya2 
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Va,2(0)  =  - 


yS 1 s 1 +YS 1  yslal  Ysls^  jS»in 

1  yalsl  yalal+y^l  yals2  ja.in 

D  ys2s 1  ys2al  ya2s2+y^2  g 

ya2al  ya2al  ya2s2  g 


This  completely  determines  the  magnitudes  Vs»^(0)  and  Va,^(0)  since  all 
the  quantities  on  the  right  side  of  Equations  (26)  and  (27)  have  been 
previously  determined. 


Note  that  the  physical  implementation  of  the  load  impedance  is 
not  important  and  in  fact  are  assumed  to  be  zero  in  the  computer  pro¬ 
gram.  They  were  included  here  in  order  to  keep  the  theory  as  general 
as  possible. 


C.  The  Transmitted  Field 

After  having  determined  the  two  voltages  Vs,^( 0)  and  Va’^(0)  of 
the  second  array,  the  transmitted  field  into  the  semi -infinite  space, 
y>d« ,  can  now  be  found.  No  radiation  takes  place  from  the  first  array 
because  of  the  shielding  effect  of  the  second  array  so  by  a  slight  change 
in  notation  of  Equation  (43)  in  [24],  the  total  transmitted  H-field  re¬ 
radiated  by  Vs,2(0)  and  Va,^(0)  is 


Y  00  00 

b  r  r 

2tnr  .  l  l 

x  z  k=-«  n=-« 


-j65(R-yb4)*r5  -jB4(yb4-R 

e _ e _ 

r5y 


-r<2). 


[in5(iP5VS,2(0)+/5Va,2(0))lT4/5(0,0) 

+  llS5(»PlvS,2(0)+,,PIVa’2(0')»T4/5(0>0)]  U8) 

where  is  the  characteristic  admittance  in  the  semi-space  y  >  b^ 
of  Figure  1  and  where  all  other  quantities  have  been  previously  defined. 

Since  the  total  transmitted  H-field  is  desired  in  the  far-field, 
all  evanescent  modes  have  disappeared.  So  only  the  k=n=0  term  i&  of 
importance  provided  no  grating  lobes  exist.  For  most  practical  appli¬ 
cations  grating  lobes  are  avoided  as  is  the  case  for  this  report.  They 
are  avoided  since  a  null  usually  occurs  near  the  frequency  at  which  the 
grating  lobe  appears. 

From  [25]  or  simply  from  conservation  of  energy,  it  is  known  that 
the  cross  polarized  component  of  the  transmitted  field  should  be  as 
small  as  possible  in  order  to  have  unit  transmission  coefficient. 
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SECTION  III 

DEVELOPMENT  OF  DATA  AND  RESULTS 


The  introduction  explained  qualtati vely  the  analysis  of  the 
structure  of  Figures  1  and  2.  Quantitative  values  are  to  be  obtained 
in  order  to  calculate  the  transmitted  field  through  the  structure 
and  hence  the  bandwidth.  To  do  so,  the  following  observations  are 
useful.  Refer  to  Figure  7  for  additional  insight. 

Assume  the  incident  field  is  in  the  YZ-plane  (a=90°).  For  an 
incident  field  orthogonal .to  the  plane  of  incidence  only  the  asymmetric 
mode  will  be  excited  (Is,in=0).  The  analysis  is  reduced  to  that 
of  an  array  of  bent  'straight'  slots  shown  in  Figure  7f.  For  an 
incident  field  paral lelatcpnthe  plane  of  incidence  only  the  symmetric 
mode  will  be  excited  (Id>  =0).  This  results  again  in  an  array  of 
bent  'straight'  slots  as  shown  in  Figure  7d.  Hence  for  at  least 
the  YZ  scan  plane  the  structure  of  Figure  1  is  similar  to  that  of 
the  structure  in  [26].  It  is  therefore  reasonable  to  use  one  of 
the  data  sets  provided  in  [27]  as  a  basis  for  the  starting  data  in 
order  to  calculate  transmission  curves.  As  the  leg  angle  L.A.  is 
decreased  from  180°  to  120°  the  above  starting  data  are  not  necessarily 
op  t i mum . 

From  [28]  it  was  demonstrated  that  a  symmetric  configuration 
for  a  monoplanar  structure  led  to  the  largest  bandwidth  of  trans¬ 
mission  curves  for  orthogonal  and  parallel  polarizations.  Using 
this  for  the  biplanar  structure  requires  that: 

d2=d4 
r'2  =e4 
ere5 

'fata  set  P27  given  in  [29]  is  then  used  as  the  initial  data  for  the 
mructure  of  Figure  1.  For  development  of  the  data  to  this  point 
iefer  to  that  reference.  Kornbau  found  that  for  an  equilateral  tri¬ 
angular  grid  (Az=D  /2)  a  leg  angle  L.A.  of  120°  yielded  the  best 
performance  with  respect  to  low  cross  polarization.  Hence  this  case 
1 s  assumed  in  the  data.  All  initial  parameters  are  now  known  and 
computer  calculations  of  the  transmission  curves  can  now  be  produced. 

It  should  be  noted  that  although  the  theory  and  computer  program 
have  been  developed  to  allow  for  any  interlace  spacing,  leg  angle 
L.A. ,  and  dielectric  constants,  they  remain  fixed  throughout.  See 
Data  (TS9)  for  values,  Table  1. 
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Since  a  direct  synthesis  is  not  feasible  at  this  time,  a  systematic 
computer  solution  approach  to  increase  the  bandwidth  of  the  transmission 
curves  was  used.  This  approach  is  to  change  one  variable  and  recalculate 
the  curves  and  keep  repeating  the  process.  This  is  possible  because  the 
effect  of  each  variable  of  the  data  set  is  qualitatively  known.  This 
allows  for  prudent  choices  concerning  changes  to  be  made  to  improve  band¬ 
width. 

The  important  variables  that  were  changed  as  part  of  the  iterative 
procedure  to  improve  bandwidth  were  the  slot  dimensions,  the  interelement 
spacing,  the  interlace  spacing  and  the  width  of  the  middle  dielectric 
layer.  The  nature  of  the  changes  follows. 

The  length  of  each  leg  of  the  slot  was  increased  in  order  to  lower 
the  resonance  frequency.  This  increased  the  bandwidth  by  moving  the  lower 
frequency  end  down  in  frequency  while  the  change  in  the  upper  end  was 
slight.  The  interelement  spacing  and  interlace  spacing  were  reduced  in 
order  to  move  the  upper  frequency  end  up  in  frequency.  The  width  of  the 
middle  dielectric  layer  determines  the  mutual  coupling  between  the  arrays. 
Since  it  was  found  by  inspection  of  Reference  (31  that  an  increase  in  the 
coupling  was  needed,  the  width  was  decreased.  Note  that  these  changes  were 
repeated  several  times. 

Each  variable  is  not  completely  independent  of  all  the  other  var¬ 
iables  so  practical  limits  do  exist.  As  an  example,  suppose  the  inter¬ 
element  spacing  is  reduced  while  holding  the  leg  length  constant.  Even¬ 
tually,  the  individual  slots  in  each  array  would  overlap.  Hence,  the 
results  would  no  longer  be  correct. 

Figures  8  to  1 1  are  the  resulting  transmission  curves  for  data  set 
TS9  using  the  voltage  distribution  given  by  Equations  (11)  and  (12). 

Notice  that  the  bandwidth  is  approaching  an  octave  for  both  orthogonal 
and  parallel  polarizations*  in  each  principle  plane  (a=0  ,90°)  for  angles 
of  incidence  up  to  75°  from  normal.  The  transmission  curves  are  in  fact 
very  similar  which  is  the  desired  result.  The  cross  polarizations,  trans¬ 
mitted  polarizations  orthogonal  to  polarizations  of  the  incident  field, 
have  been  calculated  and  are  found  to  be  -20  dB  or  less.  The  cross  polari¬ 
zations,  are  indeed  small  as  desired  to  provide  for  unity  transmission. 

After  completion  of  the  above  transmission  curves  some  develop¬ 
ments  occurred  concerning  similar  slot  elements.  Measurements  on 
a  mono-planar  slot  array  imbedded  in  an  asymmetric  dielectric  con¬ 
figuration,  A-sandwich,  showed  disagreement  with  calculated  transmission 
curves  for  frequencies  above  the  first  resonance.**  (Figures  12  and 
13.)***  It  should  be  noted  that  below  that  frequency,  the  curves  agree. 


*  This  is  in  reference  to  the  H-field  and  plane  of  incidence. 

**  Resonance  is  the  frequency  at  which  the  transmission  curve  has  no 
loss,  i .e. ,  unity  gain. 

***Figures  12  and  13  were  extracted  from  (301.  For  further  study  see 
that  reference. 
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Transmission  curves  for  orthogonal 
orthogonal  transmitted  H-field  wit 
the  voltage  distributions  of  Eauat 
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Figure  11.  Transmission  curves  for  parallel  incident  and 
parallel  transmitted  H-field  with  a=90°  using 
the  voltage  distributions  of  Equations  (11) 
and  (12).  (Data  set  TS9) 


FREQUENCY  (GHz) 

Figure  13.  Comparison  of  calculated  and  measured  transmission  curves 
for  parallel  polarization  for  a  mono-planar  slot  array 
using  the  voltage  distributions  of  Figure  12. 


^f*"»*r**.. 


The  measured  transmission  curves  did  not  contain  a  second  resonance  but 
continued  to  have  greater  loss  as  the  frequency  increased  above  reson¬ 
ance.  For  the  structure  of  Figure  1  no  measured  results  are  available. 
However,  the  same  agreement  below  resonance  and  the  same  general  dis¬ 
agreement  above  resonance  would  occur. 

To  achieve  closer  agreement,  a  new  voltage  distribution  along  the 
slot  elements  when  in  the  non-transmitting  (scattering)  mode  is  assumed. 
The  new  assumed  distribution  will  more  closely  represent  the  actual 
voltage  distribution.  According  to  [31]  the  assumed  voltage  distribu¬ 
tion  for  an  unloaded  receiving  antenna  is  of  a  cosinusoidal  form. 


VV1U) 


cosf?ds,-cosedi,ef 

l-cos^d^ef 


(29) 


instead  of  the  assumed  form  given  by  Equation  (12). 

For  frequencies  below  the  first  resonance  the  different  forms  are 
approximately  the  same.  The  transmitted  curves  below  this  frequency 
should  be  approximately  the  same.  When  using  Equation  (29)  in  a  pre¬ 
viously  developed  program  for  the  mono-planar  slot  case,  the  resulting 
calculated  curves  agreed  with  the  measured  ones  for  all  range  of  fre¬ 
quencies  (Figures  12  and  13).  This  cosinusoidal  equation  was  then  incor 
porated  into  the  computer  program  developed  for  the  structure  of  this 
report  as  the  voltage  distribution  for  the  non-transmitting  mode.  For 
the  transmitting  mode  the  voltage  distribution  is~t+re  same  as  before, 
sinusoidal.  At  this  point,  the  transmission'curves  were  recalculated. 

The  same  iterative  procedure  used  before  was  not  redone.  The  data  pre¬ 
viously  developed  was  reused  to  calculate  the  new  transmission  curves. 

The  following  Figures  14  to  17  give  the  new  computed  transmission  curves. 
It  can  be  seen  that  the  bandwidth  is  no  longer  as  constant  with  different 
incident  polarizations.  The  bandwidth  in  fact  has  been  reduced. 

A  slight  gain  also  is  observed.  We  have  not  been  able  to  explain 
the  reason  for  this  gain  although  the  transmission  curves  are  ex¬ 
pected  to  he  more  accurate  than  for  the  sinusoidal  case  above.  It 
should  be  possible  to  slightly  improve  the  bandwidth  so  that  it  is 
nearly  constant  for  the  varying  angles  of  incidence  and  different 
polarizations.  This  could  possibly  be  done  by  decreasing  the  in¬ 
terelement  spacings  and  the  interlace  spacing.  This  would  raise 
the  upper  frequency  end  of  the  transmission  curves.  Since  the  in¬ 
terelement  spacings  and  interlace  spacing  are  almost  decreased  to 
a  point  where  the  slots  overlap,  not  much  improvement  is  expected. 
However,  this  has  not  been  attempted  at  this  time. 

A  slight  descrepancy  has  arisen  when  using  Equation  (29)  as  the 
voltage  distribution  of  the  slot  elements  in  the  receiving  mode.  For 
this  passive  structure,  a  small  gain  (~  0.3  dB  or  less)  has  occurred 
for  certain  incident  angles  and  certain  polarizations.  See  Figure  16. 
This,  of  course,  is  physically  impossible. 
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and  parallel  transmitted  H-field  with  a=90°  using 
the  voltage  distributions  of  Equations  (11) 
and  (29).  (Data  set  TS9A) 


The  source  of  this  gain  has  not  yet  been  located  although  attempts 
were  made  to  resolve  the  problem.  The  computer  program  was  re-examined 
with  no  apparent  discrepancy  between  theory  and  program  found.  The 
method  of  performing  the  infinite  sums  through  use  of  the  convergence 
numbers  was  tested  without  any  changes  resulting  in  the  gain.*  Certain 
roundoff  errors  were  considered  but  to  no  avail. 

Next  the  slight  gain  was  considered  to  be  caused  by  the  method  of 
handling  the  end  effects  of  the  slots  in  the  manner  of  an  effective  length. 
Ignoring  the  end  effects  by  setting  the  effective  length  equal  to  the 
physical  length  of  each  leg,  reduced  the  gain  slightly  (gain  0.25  dB 
or  less). 

Re-examing  the  voltage  distribution  being  used  shows  that  it  was 
developed  for  infinitely  thin  slots.  The  cause  of  the  gain  was  then 
thought  to  be  caused  by  using  slots  of  finite  width.  The  width  of  the 
slots  was  then  reduced  by  an  arbitrarily  chosen  factor  of  4.  An  effective 
length  was  reinstated.  The  gain  was  reduced  (0.15  dB  or  less).  Again 
neglecting  the  end  effects  results  in  even  less  gain  (0.06  dB  or  less). 
Figures  18  to  23  show  the  critical  portions  of  the  transmission  curves 
for  each  of  the  above  cases.  Although  each  step  offered  improvement 
through  reduced  gain,  a  gain  still  existed.  Lack  of  time  has  limited 
further  search  for  the  cause  of  the  gain.  However,  the  effective  length 
does  not  seem  to  be  the  cause  since  end  effects  cannot  physically  be 
ignored.  If  the  slots  were  made  even  thinner  the  gain  would  most  likely 
reduce.  Since  the  effective  length  was  small,  the  final  transmission 
curves  were  calculated  making  use  of  it. 

Note  that  the  achieved  bandwidth  varies  from  approximately  2  GHz 
to  6  GHz  depending  on  the  angle  of  incidence  and  polarization.  This  is 
still  the  largest  and  most  constant  bandwidth  to  date.  For  angles  of 
incidence  from  normal  (1°)  to  60°  the  bandwidth  has  a  range  of  6.5  GHz 
to  10  GHz. 

Interlace  Anomaly 

Examining  the  transmission  curves  given  in  [32]  for  the  biplanar 
straight  slot  case  imbedded  in  three  dielectrics  results  in  the  follow¬ 
ing  observations.  When  the  plane  of  incidence  is  the  a=0  plane  (<fi- 
plane),  with  polarization  orthogonal  to  the  plane,  zero  transmission 
occurs  at  several  frequencies.  The  'zero'  of  interest  is  that  one  which 
limits  the  highest  frequency  that  passes  through  the  filter  with  little 
1  dB)  attenuation.  This  'zero'  will  limit  the  bandwidth  in  that  plane. 

The  cause  for  this  is  that  the  structure  anticipates  the  onset  of 
a  trapped  grating  lobe  in  the  middle  dielectric.**  When  this  occurs  the 


-3  -5 

*The  convergence  number  was  decreased  from  10  to  10  to  allow  more 
terms  in  the  summation  to  be  used. 

**Since  the  middle  dielectric  has  a  higher  permittivity,  the  grating  lobe 
will  occur  there  at  a  lower  frequency  than  in  the  other  dielectric  layers. 
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m  w=0.03  cm  and  no  end  effects 


Figure  18.  Transmission  curves  for  orthogonal  incident  and 
orthogonal  transmitted  H-field  for  the  various 
cases  used  in  discussing  the  slight  gain  that  occurs 
when  using  the  voltage  distributions  of  Equations 
(11)  and  (29).  (Data  set  TS9A  except  where  noted.) 
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Figure  19.  Transmission  curves  for  orthogonal  incident  and 
orthogonal  transmitted  H-field  for  the  various 
cases  used  in  discussing  the  slight  gain  that  occurs 
when  using  the  voltage  distributions  of  Equations 
(11)  and  (29).  (Data  set  TS9A  except  where  noted.). 
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w=0.03  cm 


Figure  20.  Transmission  curves  for  orthogonal  incident  and 
orthogonal  transmitted  H-field  for  the  various 
cases  used  in  discussing  the  slight  gain  that  occurs 
when  using  the  voltage  distributions  of  Equations 
(11)  and  (29).  (Data  set  TS9A  except  where  noted.) 


TRANSMISSION  IDB) 

-O-S  0-0  0.5  -1-0  -0.5  0.0  0-5  -1.0  -0.5  0.0  0.5 


w=0.03  cm  and  no  end  effects 


21.  Transmission  curves  for  parallel  incident  and 

parallel  transmitted  H-field  for  the  various  cases 
used  in  discussing  the  slight  gain  that  occurs 
when  using  the  voltage  distributions  of  Equations 
(11)  and  (29).  (Data  set  TS9A  except  where  noted.) 


Figure  22.  Transmission  curves  for  parallel  incident  and 

parallel  transmitted  H-field  for  the  various  cases 
used  in  discussing  the  slight  gain  that  occurs 
when  using  the  voltage  distributions  of  Equations 
(11)  and  (29).  (Data  set  TS9A  except  where  noted.) 
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w=0.12  cm 


m 


10  w=0.03  cm  and  no  end  effects 


Figure  23.  Transmission  curves  for  parallel  incident  and 

parallel  transmitted  H-field  for  the  various  cases 
used  in  discussing  the  slight  gain  that  occurs 
when  using  the  voltage  distributions  of  Equations 
(11)  and  (29).  (Data  set  TS9A  except  where  noted.) 


evanescent  waves  begin  to  increase  in  strength  and  will  destructively 
interfere  with  the  principle  wave  at  some  frequency.  This  causes  the 
mutual  coupling  betv/een  the  straight  slots  to  go  to  zero,  hence  zero 
transmission  occurs.  This  is  referred  to  as  Luebbers  anomaly  [33]. 

In  the  ct=90°  plane  (e-plane),  the  pattern  factor  limits  the  inter¬ 
ference  between  the  evanescent  waves  and  the  principle  wave,  consequ¬ 
ently  a  null  in  the  transmission  curve  corresponding  to  Luebbers  anomaly 
is  not  observed  for  orthogonal  polarization.  Luebbers  anomaly  does 
exist  in  this  plane  for  parallel  polarization. 

Since  the  starting  data  for  the  structure  of  this  report  came  out 
of  the  above  straight  slot  case,  Luebbers  anomaly  was  expected,  but  not 
observed  for  the  computed  transmission  curves.  However,  an  anomaly  ( i . e . , 
a  null  in  the  transmission  curves)  does  exist  for  parallel  polarization 
in  both  the  a=0°  and  a=90°  planes.  Investigation  of  this  anomaly  is  done 
as  before  in  the  a=90°  plane  (see  p.  19).  It  will  be  recalled  that  for 
parallel  polarization  only  the  symmetric  mode  is  excited.  (Figure  7d). 

For  a  null  in  the  transmission. curve  (Figures  15  and  17)  the  symmetric 
mutual  coupling,  Ysls2  or  Ys“'  must  go  to  zero  for  that  point  to  be  an 
anomaly.  Examining  results  of  the  computer  program  shows  this  to  be  the 
case.  In  the  a=0°  plane  both  the  symmetric  and  asymmetric.modes  are  . 
excited  but  none  of  the  modal  admittances  Ysl  ,  Y*“  ,  YalaZ  or  Ya<1 
goes  to  zero.  However,  the  anomaly  still  exists.  Hence  the  "effective 
mutual"  coupling  between  the  arrays  must  go  to  zero. 

Luebbers  anomaly  and  the  newly  attained  anomaly  are  of  different 
polarizations.  Luebbers  exists  in  the  a=0°  plane  for  orthogonal  polari- 
zation  and  in  the  a=90°  plane  for  parallel  polarization,  while  the  new 
anomaly  exists  in  both  the  a=0°  and  a=90°  planes  for  parallel  polariza¬ 
tions. 


Through  the  computer  results,  it  was  found  that  as  a  direct  result 
of  interlacing  the  slots,  Luebbers  anomaly  changed  to  the  new  anomaly. 

The  type  of  element  did  not  affect  the  anomalies.  Hence  the  new  anomaly 
is  referred  to  as  the  interlace  anomaly. 

This  interlace  anomaly  limits  the  bandwidth  at  the  upper  frequency 
ends.  As  was  previously  stated,  the  interlace  spacing  and  interelement 
spacings  were  reduced  to  raise  the  upper  frequency  end.  This  moved  the 
anomaly  up  in  frequency. 
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CHAPTER  IV 
CONCLUSIONS 


The  transmission  properties  of  a  metallic  radome  configuration 
of  two  slot  arrays  consisting  of  three-legged  elements  imbedded  in 
three  dielectric  layers  has  been  investigated  in  order  to  design  a 
bandpass  filter  with  a  large  bandwidth.  Only  the  symmetrical  case 
consisting  of  two  identical  slot  arrays  with  outer  dielectric  layers 
of  the  same  material  was  considered. 

A  mathematical  analysis  was  not  attempted  due  to  the  complexity 
involved.  A  computer  assisted  design  approach  was  used.  The  initial 
analysis  used  a  sinusoidal  voltage  distribution.  This  was  used  in  an 
iterative  process  to  obtain  as  large  a  bandwidth  as  possible.  Unfor¬ 
tunately,  this  was  found,  for  a  mono-planar  configuration,  to  produce 
erroneous  results  at  the  upper  portion  of  the  frequency  band.  Thus,  it 
became  necessary  to  introduce  a  new  cosinusoidal  voltage  distribution 
for  unloaded  three-legged  slots  in  the  non-transmitting  (scattering) 
mode.  This  should  improve  the  accuracy  of  the  transmission  curves  above 
resonance.  The  transmission  curves  for  the  design  were  calculated.  The 
bandwidth  was  then  found  to  be  somewhat  smaller  than  the  bandwidth  cal¬ 
culated  when  using  the  sinusoidal  voltage  distribution.  The  bandwidth 
was  also  found  to  vary  some  with  different  angles  of  incidence  and 
different  polarization.  For  angles  of  incidence  from  normal  (1°)  to 
60°  the  bandwidth  has  a  range  of  6.5  GHz  to  10  GHz.  Note  that  the 
iterations  were  not  repeated  when  the  cosinusoidal  voltage  mode  was 
introduced.  The  data  developed  to  that  point  was  reused.  It  is, 
therefore,  possible  that  slight  improvements  could  be  made  to  increase 
and  stabilize  the  bandwidth.  The  final  design  still  results  in  the 
largest  and  most  constant  bandwidth  to  date. 
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APPENDIX  A 

LIST  OF  ADMITTANCES  FOR  GIVEN  STRUCTURE 


In  the  main  text,  the  admittances  were  defined  and  the  general  form 
was  given.  This  appendix  lists  the  admittances  explicitly. 


yslsl  =  _  I 


slsl 


VSl(0) 


2D~r  £  £ 

2DxUz  k=-°°  n=-° 


w 

JB2  4  r2y 


2y 


tiP2Uip21iT2(°»d2)  +  llp2U/lP21||T2(0.d2)] 


2DxDz 


L  l 

k=-o°  n=-« 


■je3  T  r3y 


3y 


CiP3UiP!1iT3<0*d3)  +  ■  (A1) 

s 1  a 1  jSlal  Y2  e"jB2  4  r2y 

val(o)  2DXDZ  i^Lcc  n=-»  r2y 

tiP2UtP21iT2(0>d2)  +  „p|UHP|1(,T2(0’d2)] 

Y,  »  .  e‘jS3  Sf  r3y 

+  — - —  y  y  - 

2DxDz  k£_  nL.  r3y 

CiP3Uip31iT3((,*d3)  +  i.P3UtiP31|.T3^°‘d3^  (A2) 
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s  1  s  2  Y  a,  oo  "^3^3r3v 

Ysls2  -  _  ] _ =  _J _  i  y  e  J  J 

Vs2(0)  2DxDz  k=-<»  n=-°°  r3y 

^I’VlVs'0*0'  *  ,,P3U»PI2»T3(0-°» 


vs la2  ,  .  Is1a2  =  Y3  y  y  e 
Va2(0)  2DxDz  k=-<*  n=-°° 


■Je3d3r3y 

r3y 


t^’VfiV0-0'  +  „PlU,,P32„T3(“-0>] 


n 


(A3) 


(a4) 


To  define  the  remaining  admittances, 

change  from  to 


superscript  s 
superscript  a 


superscript  a 
superscript  s 


which  results  in  four  more  admittances. 

Then 


change  from 


to 


superscript  1 
superscript  2 
subscript  2 


superscript  2 
superscript  1 
subscript  4 


(J}T^0'd3) 


tl,T3<<l3.°> 


which  again  results  in  four  more  admittances. 

The  last  four  admittances  are  found  by  repeating  the  first  set  of 
changes.  All  admittances  for  the  structure  of  Figure  1  are  now  specified. 
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APPENDIX  B 

T-FACTOR  FOR  THE  STRUCTURE  OF  FIGURE  1 


Note  again  that  a  slot  can  be  considered  as  a  magnetic  element  in 
front  of  an  electrically  perfectly-conducting  ground  pTane.  The  re¬ 
flection  coefficient  for  tne  H-field  from  a  ground  plane  is  one.  When  cal¬ 
culating  self  admittances,  the  remaining  arrays  are  short  circuited 
causing  them  to  become  simply  ground  planes.  Therefore,  for  either  self 
or  mutual  admittance  calculations  the  reflection  coefficient  at  each  slot 
array  become^  one.  F^r  the  structure  of  Figure  1,  the  reflection  coeffi¬ 
cients,  .123,  x  1  32*  fi,l34  and  j^3  are  equal  to  one. 

For  the  self  admittance  calculations  the  reference  point  of 
the  single  test  element  is  considered  to  be  one  quarter  the  slot  width 
away  from  the  reference  point  of  the  array  by  definition.  However, 
the  T-factor  is  calculated  assuming  the  reference  element  and  the  array 
coincide.  This  discrepancy  is  negligible. 


For  mutual  admittance  calculations  the  reference  element  coincides 
with  one  of  the  arrays. 

Using  the  above  observations  in  the  generalized  non-normal i zed  T- 
factor  (Equation  (17))  results  in  the  following  for  the  given  structure 
for  self  admittance  calculations; 


±  T2(0,d2)  =  2.0 


v  'Hi1*-' 


-j2f52d2r2y\ 

e 

^-j2e2d2r2y 


(Bl) 


T3(0,d3)  =  2.0 

lli; 


L  +e-^3d3r3y\ 
Vi  .  e'j2B3d3r3yJ 


(B2) 


i  T3(d3«0)  =  .  Tt( 0,d3) 

f{>  3  (h!  3  3 


(B3) 
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and  for  mutual  admittance  calculations 


{i}T3(0,0)  =  4/(l  - 


■j^3d3r3y 


■) 


( 


This  provides  all  necessary  T-factors  for  admittance  calculations. 
The  normalized  T-factors  needed  for  the  structure  of  Figure  1 


,mt(0,0)2/,  =  2 


1  ~  (ji)ia.i _ 

~  '"’p  -J^2d2r2y 


and  for  transmitted  H-field  calculations 


{i}T(0,0)4/5  =  2 


1  ’  T 


p  -J'234d4r4y 
1  '  {A}  4,5  e 


ATI  T-factors  have  now  been  determined. 


APPENDIX  C 

THE  PLANE  OF  INCIDENCE  AND  THE  PLANES  OF  SCATTERING 


The  following  formulas  were  taken  from  Appendix  B  in  [34]  and 
are  generalized  and  included  here  for  completeness. 

The  plane  of  scattering  is  defined  as^the  plane  containing  the 
vector  normal ^to  the  dielectric  interface,  nQ,  and  the  direction  of 
propagation,  rm,  where  m  refers  to  the  dielectric  media. 

For  unit  vectors,  inm,  orthogonal  to  the  plane  of  incidence 


inm 


noxrm  _ 

IWnl 


■xrmz+zrmx 


(r  z+r 
v  mx  mz  ' 


2x1/2 


(Cl) 


For  unit  vectors  parallel ,..0^,  to^the  plane  of  scattering  and  orthogonal 
to  the  direction  of  propagation,  rm. 


n  nm  =  inmxrm 


1 


(rmv2+r_/) 
'  mx  mz  ' 


s-Ti,  (-xr  r  +y(r  2+r  2)  -  zr  r  ) 
2^1/2  mx  my  mx  mz  my  mz 


(C2) 


Thg  plane  of  incidence  is  defined  as  the  plane  containing  the 
vector  nQ  normal  to  the  dielectric  interface  and  the  direction  of 
propagation  s  (=r  for  k=n=0),  cf.  -F-igure  Dl. 
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APPENDIX  D 

REFLECTION  COEFFICIENTS  FOR  THE  H-FIELD 


The  following  formulas  are  found  in  Appendix  C  in  [35]  and  are 
listed  here  for  completeness.  The  formulas  correspond  to  the  following 
generalized  figure. 


Figure  D1 .  Incident  wave  on  a  dielectric  boundary. 


fz 


m+1 


(m+1 )  rmy  *  ’Sn1"  (ro+i  )y 


/ 


e  (m+1 )  rmy+>/^mr  (m+1  )y 


11  Hn  ,m+l 


/ 


e(m+l)r(m+l)y~/Emrmy 
^e(m+l )r(m+l )y+/^mrmy 


l+jJ""m,m+l  ^c(m+l  )rmy 


(Dl) 


(02) 


(03) 


mr(m+l )y 


APPENDIX  E 

MUTUAL  ADMITTANCE  (Yas  AND  Ysa)  BETWEEN  SYMMETRIC 
AND  ASYMMETRIC  MODES 


In  this  appendix  we  assume  that  the  voltage  distribution  on 
the  slots  is  of  the  same  form  for  both  the  transmitting  and  non-trans¬ 
mitting  modes,  namely  sinusoidal.  The  admittance  relationship  presented 
also  holds  for  the  cosinusoidal  voltage  distribution  around  resonance 
where  the  two  voltage  distributions  are  approximately  equal.  Any 
interlace  spacing  is  allowed,  i.e.,  the  relationships  hold  for  straight 
or  interlaced  array  grids. 

The  relationships  that  follow  are  developed  so  that  easy  compari¬ 
son  of  computer  program  results  can  be  made  between  Ysa  and  Yas. 
Separating  Ysa  into  real  and  imaginary  parts  due  to  real  and  imaginary 
space*  yields 


/sa 


=  (ReYsa+ImYsa)real  +  (ReYsa+ImYsa) 
space 


imaginary 

space 


(El) 


Similarily 


Yas  -  (ReYas+ImYas)  ,  +  (ReYas+ImYas) . 

real  imaginary 


(E2) 


space 


space 


Equating  the  components  to  each  along  with  extensive  use  of 
Equations  (B5),  (B7)  and  (B8)  in  [36]  it  was  found  that 


ReY 


sa 


imaginary 

space 


ReY 


as 


imaginary 

space 


(E3) 


ImYsa.  .  sa 

imaginary  ~  +  I™/  imaginary 
space  space 


(E4) 


♦Real  space  pertains  to  propagating  mode,  i.e.,  k=n=0  for  no  grating 
lobes.  Imaginary  space  pertains  to  evanescent  modes,  i.e.,  k  and/or 
n/0  for  no  grating  lobes. 
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for  a  three-legged  element,  provided  that  the  T-factor  is  real  for 
imaginary  space.  This  holds  if  the  effective  reflection  coefficient 
is  real,  namely  if 


1)  no  grating  lobes  exist,  anywhere,  or 

2)  no  dielectric  exists  (i.e.,  free  space). 

An  additional  case  occurs  if  the  reflection  coefficient  is  complex 
but  multiplied  by  a  negligible  number,  effectively  making  the  T-factor 
real.  This  is  the  case  if 

1)  the  electrical  thickness  B_.d_.r-  ,  of  the  dielectric  slabs 
is  sufficient  to  insure  thSt^tFre  exponentials  in  the  T- 
factor  be  negligible. 

Equations  (E3)  and  (E4)  correspond  to  Equation  (B18)  in  [37]. 

For  the  three-legged  elements 


)  =  z 

(E5) 

p(2)  =  i  p(2)+£  p(2) 

(E6) 

p(3)  =  -  x  p(2)+z  p(2) 

r  V  7 

(E7) 

the  relationship  between  Ysa  and  Yas  in  real  space  is 

(Ysa)real  *  ‘  <vaS>real 
space  space 


or 


★Note  Equation~(E8)  and  Equation  (B16)  in  [38]  have  differently  defined 
quantities.  Hence  they  do  not  contradict. 
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<Rev“)rea,  *  (taYsa)rea1  -  -  (Rev^)rea1  -  umvas)real  •  (E9) 

space  space  space  space 


For  the  above  relationship  to  be  true  it  was  found  that  for  the  sym¬ 
metric  pattern  factor,  Ps,  the  x  and  z  components  had  to  be  pure 
imaginary  and  pure  real,  respectively,  and  for  the  assymmetric  pattern 
factor,  Pa,  the  x  and  z  components  had  to  be  pure  imaginary  and  pure 
real  respectively.  This  must  also  hold  for  the  transmitting  pattern 
factors.  For  the  given  three-legged  elements  this  is  indeed  true  only 
when  the  incident  field  is  in  the  principle  planes  (a=0°,90°). 


51 


APPENDIX  F 

PATTERN  FACTORS  IN  THE  YZ-PLANE 
FOR  INTERLACE  STRUCTURE 


This  appendix  assumes  that  the  voltage  distribution  on  the  slots 
is  of  the  same  form  for  both  transmitting  and  non-transmitting  modes. 
Certain  properties  of  the  composite  pattern  factors, ri;Pst,  riiP  , 

1  ii 1  Mr 

,,,PS  and  n,Pa,  will  be  established  for  interlacing  in  the  z  direction 
Mr  Mr 

only  (ax=0).  The  properties  of  the  composite  pattern  factors  for  the 
non-interlace  structure  are  a  special  case  of  those  for  the  interlace 
structure  treated  here,  (see  Appendix  D  in  [39]. 

When  the  plane  of  incidence  is  the  YZ-Plane,  smx=0 


--  +  y  r  +  z 
Dx  J  my 


(FI) 


where  subscript  m  refers  to  the  media. 

Using  a  three-legged  element  that  is  symmetric  with  respect  to  the  YZ 
Plane,  results  in 


pO)  =  z  (F2) 

p(2)  _  x  p(2)  +  z  p(2)  (F3) 

p(3)  =  -  x  p^2^  +  z  p^2^  (F4) 

Since  the  overall  objective  is  to  find  certain  properties  between 
the  composite  patterns,  certain  properties  between  the  patterns  of  each 
leg  need  be  found. 

To  find  those  properties  the  following  Equation  (F5)  is  assumed 
since  this  relates  the  exponential  factor  of  the  patterns  of  each  leg. 
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1 


P^-rm(k2,n)  =  p(3^rm(k3,n)  (F5) 

where  k^  and  k3  refer  to  the  summation  indice  k  for  leg  2  and  3, 
respectively. 

Substituting  Equation  (FI)  into  Equation  (F5)  gives 


so 


(F7) 


or 


2nAz 

Dz 


k  -  k 
3  2 


(F8) 


The  relationship  between  ( k£ » n )  and  rmx(k3»n)»  where  rmx  is 

the  x  component  of  rm,  can  now  be  found.  '' 


or 


rmx^2’n^  "  rmx^3,n^ 


(F10) 


Also  by  inspection  r^^.n^r^^.n). 

Substitution  of  Equations  (F5)  and  (F8)  into  the  equations  for 
the  pattern  results  in  an  equation  of  the  form 
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pv2  ^2rmz  _  k3 ,n)  =  Pv3(k3,n) 


(FIT ) 

Pv3  ^  .  k3>nj  =  Pv2{k3>n) 

where  Pv  is  the  symmetric,  asymmetric,  non-transmitting  (scattering) 
and  transmitting  patterns,  Ps ,  P3 ,  Ps*,  P3*  of  each  leg  for  scan  in 
the  YZ-Plane.  v  is  a  dummy  variable. 

Substituting  Equation  ( FI  1 )  into  (DIO)  in  [40]  with  slight  notational 
change  yields 


lPm(k2’n)  9  o  1/2 

(r2x(k2,n)+r2z(k2,n))1/2 

[2rmx{k2.n)pSl(k2)n)+px>"mZ(k2*n)(pS2(k2»n)"pS3(k2’n^ 

-pz  rmx(k2»n)(pS2(k2»n)+Ps3(k2,n)] 


1 


<rL<k3'n)trmz<k3-n>>’/2 


-2''m*(k3-")pSl(k3-">*Pxrmz<k3-n>  -k3.")  -^(tT  -k3’")) 


+  p  r 
z  mx 


(k3,n)  (ps2(2jjaz  .^.n)  *Ps30nJZ  -k3,„)) 


(FI  2) 


1 


<rmx<kf>n)+rmz<k3"n>)1/Z 

t-2rmx(k3,n)pSl(k3in)+pxrinz(k3>")(pS2(k3,,')*pS3,k3>") 

+pzrmx ( k3  »n) ( pS^ ( k3 ,n)+pS3 ( k3 ,n) ) ] 


"  -  lrk3’^ 
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Substituting  Equation  (F7)  into  Equation  ( FI 2 )  yields 

ip„  (Tf  -k3'n)  ’  ipm<k3-">  (F'3) 

Similar  derivations  using  Equations  (Dll )-(D13)  in  [41  J  of  non-trans¬ 
mitting  modes  for  parallel  and  orthogonal  components  yields 

"p™  (Tf  ~k3’n)  *  »pm(k3'")  <F'4> 

<n5> 

»Pm  ~k3’n)  ■  -  „pm<k3'">  <F'6) 

The  above  for  equations  also  hold  for  the  transmitting  case. 

Note  that  for  k=n=0. 

^(0,0)  =  nPa(0,0)  =  0 


which  is  identical  to  Equation  (D18)  in  [42].  This  should  be  the  case 
since  at  k=n=0  the  principle  propagating  mode  will  be  the  same  regard¬ 
less  of  grid  structure. 

If  az=0  (non-interlace  structure)  Equations  (F13),  (F14),  (F15) 
and  (F16)  reduce  to  Equations  (D14)-(D17)  in  [43]  given  for  the  non¬ 
interlace  design. 

Equations  (D21)  and  (D22)  in  [44]  only  hold  for  the  non-interlace 
structure. 
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APPENDIX  G 

COMPUTER  LISTING  FOR  BIPLANAR  SLOT  ARRAY 
OF  THREE-LEGGED  ELEMENTS  IN  A 
STRATIFIED  DIELECTRIC  MEDIUM 


LIST  OF  IMPORTANT  EXPLANATION  AT  PROGRAM* 

COMPUTER  VARIABLES  LINE  NUMBER 


ABRS .  91 

AC  PAT .  658 

ALPHA .  17  DATA  FILE 

APATX .  641 

APATZ .  642 

BETA .  85 

BETAD .  891 

BTD .  586 

CA1 .  1202 

CS1 .  1201 

D .  11  DATA  FILE 

DELTZ .  26  DATA  FILE 

DX .  24  DATA  FILE 

DZ .  25  DATA  FILE 

EFFL .  893 

ELEMX .  894 

ELEMZ .  894 

ER .  19  DATA  FILE 

ETA .  17  DATA  FILE 

EXPT .  1127 

EXPYM .  1103 

EXPYS .  1077 

FINCRM .  9  DATA  FILE 

FREQH .  8  DATA  FILE 

FREQL .  7  DATA  FILE 

IKK .  853 

INN .  853 

IOC .  52  PRINT 

IVD .  888 


*  Except  when  noted  otherwise.  Explanation  then  exists  at  line 
number  given  in  listed  program. 
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LIST  OF  IMPORTANT  EXPLANATION  AT  PROGRAM 

COMPUTER  VARIABLES  LINE  NUMBER 

(cont.)  (cont.) 


NADMT .  587 

NANGLS .  15  DATA  FILE 

OHO .  1221 

OHP .  1222 

OPATC .  660 

PATX .  897 

PATZ .  897 

PHO .  1223 

PHP .  1224 

PPATC .  661 

RHO .  1042 

RL .  27  DATA  FILE 

RLA .  30  DATA  FILE 

RLAMDA .  85 

RLF .  585 

RXKN .  859 

RYKN .  859 

RZKN .  859 

SCPAT .  624 

SPATX .  637 

SPATZ .  640 

SX .  105 

SZ .  105 

TFACT .  997 

THIK .  29  DATA  FILE 

VA2 .  1251 

VS2 .  1250 

WIDTH .  28  DATA  FILE 

YADMT .  590 

YADMTL .  591 

YADMTR .  592 
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DATA  FILE 


Figure  Gl.  Program  structure. 
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Figure  G2.  Defining  the  physical  structure. 
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